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INTRODUCTION 


A great deal of interest has been generated recently in the possibility of producing new materials 
hi the reduced gravity environment provided during the forthcoming missions of Spacelab. The range 
of possibilities extend from producing large crystals of uniform properties to manufacturing materials 
with unique properties. Most of these processes involve the solidification of materials from the liquid 
state. Convective motions within the liquid during solidification can influence the local material 
composite and the shape of the solid-liquid interface which may result in solids with non-uniform 
properties and crystal defects. The microgravity environment of Spacelab is being viewed as one in 
which the buoyancy forces are eliminated so that convection driven by thermal gradients will not 
occur. It is hoped that this will lead to an improved solidification process. However, convection may 
occur for other reasons and whether convection is negligible or not during solidification constitutes 
processing in low-gravity environment. Little information exists presently on convection during 
solidification under such circumstances. 

The work reported here is a continuation of an analytical investigation into the nature of 
convective motion in a binary liquid layer due to surface tension forces during its solidification. The 
onset of convection will be determined through > ->• ability analysis which is described below. 

STATEMENT OF THE PROBLEM 

The occurrence of convective motions in a Siu i has been studied both theoretically and ex- 
perimentally for approximately a century. The problem is very well documented in several books 
and numerous articles with all of its variations. It is obvious that in order to study analytically 
the detailed convective motions in a fluid in any configuration requires nothing less than the total 
solution of the Navier-Stokes equations and Energy Conservation equations. This is well known to 
be a formidable, although not an impossible, task requiring considerable amounts of time as well as 
financial resources. 

In certain applications, it is sufficient to be able to know whether or not under certain conditions 
a fluid could sustain convective motions. The answer to such a question requires far less work; 
but, again, the information is essentially just of a binary form. Such information can be obtained 
through hydrodynamic stability analysis. Essentially one introduces perturbation into a well-known 
basic state and studies the evolutions of these perturbations in time. If the perturbations decay with 
time, then the basic state is said to be stable and no convective motion will ensure. If, on the other 
hand, the perturbations are found to grow with time, then the basic state is said to be unstable 
and convective motion will take place. Fundamental to such a stability analysis is the existence of a 
stationary basic state. Unfortunately, the problem under consideration does not have such a basic. 
However, it has been shown that it is possible to carry out a meaningful! stability analysis of such 
a basicc state through several techniques. We have chosen the energy stability method as the best 
suited for the problem under consideration. This technique is elucidated in detail in the monographs 
of Joseph [1]. 

In this work, we consider the stability of a binary fluid layer which is being solidified from below 
and has an upper free surface. The fluid layer is assumed to be of infinite extent in the horizontal 
direction. Since it is assumed that the process is being carried out in the low gravity environment of 
Spacelab, it is anticipated that the driving force in the surface tension force at the free surface. The 
solution to the basic state has already been obtained in a previous report (Antar [2]). As is expected, 
the basic state is a function of both time aud space. Thus, the stability analysis used must account 
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for the variations of the baste state with time. 

The stability analysis starts with the perturbation equations which may be written in the 
following form: 


+ + + + fori\k (l) 

at p 0 

+ .V7 + tt- VC = DV 2 7 (2) 


+ + tt.VT = /TV 2 * (3) 

where u (n, y, w) is the perturbation velocity vector, U (U, V, W) is the basic state velocity vector. 
7 and C are the perturbation and basic state concentration and 9 and T are the perturbation and 
basic state temperature, respectively, u, and D and k are the kinematic viscosity, thermal diffusivity 
and the solutal diffusion coefficient, respectively. These equations are for a Boussinesque fluid and 
subject to the following boundary condition 


« = 0 = 7 = 0 


on the lower surface. 

While at the upper surface, we must have the following boundary conditions: 


.at 

(4a) 

£-0 

dz 

(46) 

t dw . du . da 

' ,( 8? + 9l )= 5I 

(4c) 

,dw . dv . 9a 

m 


where a is the surface tension. On all the sidewalls we must have the perturbation function vanishing. 

The governing equations are first non-dimensionalized using the fluid depth d as a length scale, 
dPfn for a time scale, n/d for a velocity scale, T\ — T 0 and Ci—C t for temperat ure and concentration 
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scales respectively. Upon introducing these scales into equations (1) - (4) we get the following 
equations for the perturbations functions: 


~ + u • Vti + U • Vu + u • VU « - Vp + PV 9 u + P(R6 + LRci) 

ut 


(5) 


M 

dt 


+ U-V0 + U‘V0 + u-VT=* V a 0 


^7 + U- V7 + u. V7+«.VC = LV 2 7 
at 


( 6 ) 

(7) 


with the following boundary conditions: 


0 = u = 7 = O (8) 

on all side except the upper surface and 


|j + S» = 0 

(9a) 

o 

II 

|<S 

m 

du . . ,d$ . .O') . 

_ + m_+^l 5 1 = o 

(9c) 

^- + M^- + M c Lp-^0 

oz dy dy 

m 


at the upper surface. In these equations we have the following non-dimensional numbers. 
P = Prandtl No. = £ 

R Rayleigh No. = fiigST^ 

R e = Solutal Rayleigh No. = gPtfC-jfij 
L = Lewis No. = £ 

B = Biot modulus = ^ 

M = Marangoni No. (=ffi)6TdnK 
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M c — Solut&l Morangoni No. **= 

■where jf and |£ are the constants from the surface tension variations -with both temperature and 
concentrations. 

Now to obtain the energy equations we first take the dot product equation (5) with u and 
integrate about the volume under consideration to get 

^ - f P~ l q?/2dv a* — f V# : Vudv + R f w$dv -f LR t f w^dv 
dt Jv A *v Jv 

—M f 0—2-dxdy — M e L f 'l^-dxdy (10) 

■l|B 1 OZ J*<rn 1 02 

Similiarly it is possible to obtain equations for 0 2 /2 and r i 2 /2 in the following integral form 


4 - f Pfidv = - / 1 u$^dv -B [ fidxdy ~ [v$- VOdv (11) 

dt J v Jff OZ Jm~1 Jv 


it /, — /. - l /. ^ ^ 


( 12 ) 


Now multiplying equation (11) by X# and equation (12) by X 7 and adding equations (10) - (12) we 
fet 


j t K « -D + J[(R - + (R e L - \^~)wy]dv 

- f \M9 + M c Ll$rd*dy (13) 

Jgaml OZ 

where 

K - | J(P ~ + X,* 2 + 

D = I [Vu : Vu + \,VB • V6 + X^LV^ • Vi]dv + B f \ t Pdxdy 
Jv *»™1 

Equations (13) can be cast in a symmetric form with the following change of variables 


f = y/\fO\ X' = y/\^ 
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\$ ass Mfi$; X, sat Men*! 


R - MN;R e = M e N c 

'which after dropping the primes takes the following forms: 


wkwe 


and 


dflT 


-P + /MII 


K- l -f(P-'i‘ + P + '?)4v 

Z?«= /(Vo : Vu + V0- V0 + V^- V7)dv + B / 

** *S— ] 


n as CXtfiY 2 ( j CSwddv — J 


h dw 


C6w8dv — / 0~didy) 
oz 


+ct’,nl/*{f v C' t vndv - n^dxdy), 


Q$ 


ft 


[Me 

WaT 


C^LN'-n, 


dT 


6z 

9C 

dz 


0 2 dxdy 


M =M + Mc 
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Now let us focus our attention on Eq. (14). Since then, the energy of the perturbation will 
increase or decrease with time depending on whether y/kUfD^rO. Of course the direct and only 
way of gaining this information is by solving the integro-differential system explicitly. However Eq. 
(14) can be written as an inequality in the following form: 


1 dK 
D it 


-l + v/>f|:< -1 + 


SM 

P 


(15) 


where 


p * 1 »Fn 


(16a) 


with 


0 = 1 


(16ft) 


where now the problem is cast as a variational problem defined by (15) and (16). 

Let us introduce the lagrange multipliers M^t) and p(x, y, z, t ) into the variational problem (15) 
• (16) with the following optimization procedure resulting 

a(n+ /^ (v '“ )dt ' - ]5; D)=0 <n) 

From variational calculus we know that the solution to (17) is the same as the solution of the 
Euler-Lagrange equation resulting from (17). For the present case, the Euler-Lagrange system of 
equations are the following 


~ v 2 « = o 
dx 

(18) 

|E ~ v 2 v = o 

dy 

(1®) 

Qt M *c,$+ M 'C n r 9 a P + V 2 to = 0 

y/U, • '/Ml 1 d ‘ ■ 

(20) 

Oh. /S£rC,w + r s » — o 

Z yfiL, 

(20 
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^!LJ!2-C 7 w+V a T- 0 
* Vm, 7 


( 22 ) 


8u , do 
9t + 9v 



(23) 


with the following boundary conditions 


8z 2 m 


at 89 
y/i*t dx 



o-r 8n 
y/ii^dx 


= 0 


(24o) 


8v_ , l M at Of 

dz' T 2 l *y/]i t dy 


+\ M > 


a 7 fly 


0 


(246) 


88 . 1 . . o# 8w 


1*8 + % + ^ 


yfii'dz 


(24c) 


07 1 i 1/ 6*1 Ow 

al + 


0 


(24d) 


and z = 1 and 




(24c) 


at z = 0. 

Now as discussed earlier since the basic state itself a function of time then the term II in the 
energy equation (14) is a function of time. Also, the lagrange multipiers and the term in (16a) are 
also function of time. Thus, the variational problems (15) - (16) must be solved at each instant of 
time which requires the solution of the Euler-Lagrange equations at each point in time. Thus, the 
procedure is to solve the system (18) - (24) at each instant of time and obtaining the optimization 
parameter (t) as a function of time. 
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